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Abstract - We study the flow of bright solitons through two asymmetric potential wells. The 
scattering of a soliton by certain type of single potential wells, e.g., Gaussian or Rosen-Morse, is 
distinguished by a critical velocity above which solitons can transmit almost completely and below 
which solitons can reflect nearly perfectly. For two such wells in series with certain parameter 
combinations, we find that there is an appreciable velocity range for which solitons can propagate 
in one direction only. Our study shows that this directional propagation or diode behavior is due 
to a combined effect of the sharp transition in the transport coefficients at the critical velocity 
and a slight reduction in the center-of-mass speed of the soliton while it travels across a potential 
well. 
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' ■ Introduction. — The interest in the problem of soli- 
^-Htons scattering by external potentials such as barriers 
; pTpIKZ], wells [SHU], ste P s [SHU], and surfaces [HUH], 
i stems from the nonlinearity involved and the wave nature 
' C*~)of solitons characterizing such process. For instance, the 
. nonlinearity in the Schrodinger equation that controls the 
' evolution of solitons, leads to the nonclassical behavior of 
• • partial transmittance and partial reflectance upon scatter- 
' .^ing a soliton by a potential barrier. The wave nature of 
solitons results also in a host of other nonclassical scat- 
?— i tcring outcomes. Of particular importance is the so-called 
i . 4 uan tum reflection [5D] . It has been shown that for slow 
enough solitons, complete reflection from a potential well 
|13| or a down potential step takes place |17j . Quantum 
reflection is caused by a transport of matter or energy from 
the incident soliton to the bound state of the well which 
interacts repulsively (attractively) with the original soli- 
ton if these are out of phase (in phase). Furthermore, it 
was found that for a certain type of single potential well, 
e.g., Gaussian or Rosen- Morse (RM), the two scattering 
regimes, namely complete quantum reflectance and com- 
plete transmittance, are separated by a sharp transition 
taking place at a critical velocity [13] , and this interesting 
phenomenon has recently been explained for other systems 
using a delta potential [211122) . 

Another important effect is the reduction in the center- 



of-mass speed of the soliton as it crosses the potential re- 
gion. This is due to the fact that transmittance is not 
complete; a small fraction of the soliton is reflected while 
most of the soliton is transmitted. The reflected part is 
typically dispersive and decays in the form of radiation. 
The larger fraction, i.e., the transmitted part, retains its 
soliton shape and propagates but with some collective os- 
cillations. The energy for this kind of internal excitation 
comes from the center-of-mass kinetic energy. As a result, 
the center-of-mass motion of the soliton slows down. 

These two effects, the criticality and the velocity reduc- 
tion, can be combined in a useful way when two or more 
RM potentials are placed in series and can lead to novel 
phenomena. In this letter, we study the dynamics of a 
soliton in the presence of two RM potential wells with 
slightly different parameters. One of the interesting ef- 
fects of this potential combination is the directional flow 
of solitons. For certain parameter combinations, there is 
a velocity window for which solitons can propagate freely 
with nearly full transmission in one direction and blocked 
with almost perfect reflection in the other direction. We 
termed this behavior as "soliton diode" . Asymmetric elec- 
tron flow with two tunneling barriers was studied [53 26 
and for solitons this behavior has recently been demon- 
strated for a discrete nonlinear Schrodinger equation with 
asymmetric potentials [27] . While we obtain the transport 
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coefficients for the soliton diode numerically by solving the 
appropriate Gross-Pitaevskii (GP) equation, we support 
our qualitative explanation of this effect by investigating 
in detail the dependence of the critical velocity and the ve- 
locity reduction for each potential well separately in terms 
of the parameters of the wells and the speed of the inci- 
dent soliton. This allowed for an independent and accu- 
rate account of the velocity window for which this diode 
functions. 

Theoretical Model. — In the presence of an external 
potential V(x) the dynamics of bright soliton is governed 
by the one-dimensional GP equation, which can be writ- 
ten in standard dimensionless (for explicit units, see, for 
example, Ref. [55]) form as 



-\^ + V{ X )+ 9 \^{xM 



tj}{x,t), (1) 



where we choose the potential V(x) of the form 

V(x) = —V\ sech 2 [ai(x - x cm i)]—V 2 sech 2 [a 2 {x - x cm2 )]. 

(2) 

Here, Vi f2 , «i,2, and x C mi,2 determine the depth, inverse 
width, and the position of the center of the first and sec- 
ond potential wells, respectively. This potential, which 
is known in the literature as the Rosen-Morse potential, 
belongs to the class of rcflcctionless potentials which may 
transmit without reflection in the linear regime. The fac- 
tor g denotes the mean-field interaction strength and we 
take g = — 1. In our numerical simulation we prepare an 
initial state far away from the potential region and prop- 
agate it in real time. We always choose our initial wave 
function as the exact solution of the homogeneous version 
of Eq. |T]), namely 

il>(x,t = 0) = Ae iv ( x - x ^sech[A( X -x cm )}, (3) 

with A = 1 and x cm = xq + v t, where v and xq are the 
initial center-of-mass velocity and position, respectively. 
The evolution time is taken long enough for the scattered 
soliton to be far away from the potential. Then, we calcu- 
late the reflection, trapping, and transmission coefficients. 
For the left-to-right moving soliton with a single poten- 
tial well, we define the reflection, trapping, and the trans- 
mission coefficients as: R = (l/N) J _ \ip(x, t)\ 2 dx, L = 

f_ t \ip{x,t) | 2 dx, and T = (1/N) |^(x, t)\ 2 dx, 
where I sa 5/ a from the center of the well and N is the nor- 
malization of the soliton given by N — J_ \ip(x, t)\ 2 dx. 
For right-to-left moving soliton, R and T are interchanged 
but L remains the same. For the present case of two po- 
tential wells in series, we take —I to be on the left of the 
left well and I to be on the right of the right well. 

Results: The Diode Behavior. — We present the 
main findings of our study here. We consider two RM 
potential wells in series with slightly different parameters 
(See the schematic figures in Figs. [Ha) and (b) below.). 



For this set of potential wells, the soliton is reflected for 
small enough velocity and transmitted for large velocity. 
The transition from high reflection to high transmission 
is very sharp. Figures [T] (a) and (b) show the reflectance 
R and transmittance T versus the magnitude of the inci- 
dent soliton velocity v for the potential given by ^ with 
V\ = 4.35, V 2 = 4.0, a\ — a 2 — 2, and x cm i = -x cm2 = 6. 
Solid and dashed lines in Fig. QJa) show the reflection co- 
efficient for the right-to-left moving and left-to-right soli- 
ton, respectively. Figure [T](b) shows the transmittance for 
both cases. For the left-moving soliton there is a critical 
velocity for which there is a sudden jump from mostly re- 
flectance to mostly transmittance. Similar behavior exists 
for the right-moving soliton with another critical velocity. 
Surprisingly, there is a difference in between the two crit- 
ical velocities for the left- and right-moving solitons. As 
a result, there is a small but appreciable velocity window 
for which there is almost full transmittance in one direc- 
tion and nearly zero transmittance in the other direction, 
i.e., the soliton shows directional propagation for this set 
of parameters. This behavior is similar to the diode effect 
in semiconductor physics. Figures OJa) and (b) show the 
spatio-temporal plots corresponding to a selected value of 
the initial speed. FigureUJa) corresponds to a soliton inci- 
dent from the right of the right well with incident velocity 
v = —0.33. The soliton is first transmitted through the 
right well (centered at 6) and then reflected to the right 
by the left well (centered at —6) and finally transmitted 
to the right through the right well. The overall reflection 
coefficient is R 0.98 and is marked by a filled circle in 
Fig. [lja). Figure [Hb) corresponds to a soliton incident 
from the left of the left well and propagates to the right 
with the same magnitude of the incident velocity v = 0.33. 
This time the soliton is transmitted through both wells. 
The overall transmission coefficient is T ps 0.99 and is also 
marked by a filled circle in Fig.QJb). 

This novel diode behavior of soliton propagation can 
simply be explained by a combined effect of the critical 
velocity and the reduction in the soliton's center-of-mass 
speed caused by a potential well. Consider a soliton in- 
jected on a single potential well with a speed slightly larger 
than the critical velocity. The flow is naturally symmetric 
with full transmission from both sides. Introducing an- 
other potential well with considerably lower critical veloc- 
ity, say to the right of the first one, breaks this symmetry. 
While the soliton incident from the left transmits both 
the wells, the soliton from the right with the same inci- 
dent speed can only transmit the well on the right since, 
as it crosses that well, it suffers a velocity reduction which 
is small but enough to reduce the soliton's speed to below 
the critical velocity of the left well. Therefore, the soliton 
will be reflected from the left well to the right. Since the 
critical velocity of the right well is still considerably lower 
than the speed of the reflected soliton, the soliton will 
continue its motion in the right direction. In this man- 
ner, both the initially left-moving and the right-moving 
solitons end up moving to the right. 
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As a specific example, we calculate, for the potential 
wells considered in Figs.^a) and (b), the critical velocities 
of the left and right wells, v cr and v^ r , respectively, which 
take the values: v l cr = 0.3275 and v r cr = 0.1175. Let 
us consider the case of the soliton with incident velocity 
v = 0.33. For this incident velocity, the velocity reductions 
are Av l = 0.015 and Av r = 0.014. If it is incident from the 
left, it will transmit through the left well as the velocity 
is larger than v cr . The transmitted soliton will have the 
reduced velocity v = 0.315 which is still larger than v r crl 
allowing the soliton to transmit through the right well to 
the right. Now if the soliton is incident from the right side 
of the right well, it will transmit through it, as the soliton's 
velocity v = 0.33 is larger than v r cr . Due to the velocity 
reduction by the right well, the soliton's velocity reduces 
to v = 0.316 which is less than the v l cr , and therefore, 
the soliton will be reflected by the left well to the right. 
Finally it will be transmitted through the right well to the 
right. So whether the soliton is incident from the right or 
from the left, it will always end up moving to the right 
and there is no flow in the left direction for the velocity 
considered. 

Having demonstrated the function of the diode for a spe- 
cific set of the system's parameters and a single soliton's 
initial speed, a thorough investigation of this phenomenon 
in terms of these parameters is in order. Specifically, we 
want to account for the velocity window in Fig. [T] over 
which the diode functions. In addition, we want to find 
the effect of the potentials' parameters on the width of this 
window with the aim of finding the optimized potentials' 
parameters that lead to the maximum velocity window. 

We start by characterizing the effect of a single potential 
well, V(x) = — Vo sech 2 (aa;), on the propagation of soli- 
tons. We found that it is essential to set a = a/Vo' in order 
to get a sharp transition of the transport coefficients at 
the critical velocity. For such particular type of potential 
well, the width of the well decreases with increasing the 
well depth. The critical velocity turns out to be a mono- 
tonic decreasing function of the well depth. Figure [2ta) 
shows the critical velocity, defined by v cr = (vi + v 2 )/2 
as a function of the well depth Vo- Here v\ is the largest 
velocity with almost full reflection and v 2 is the smallest 
velocity with nearly perfect transmission. Circles show 
our numerical data and the solid line shows the fit of the 
data by v cr ps 0.42 V" " 018 . Currently we are unable to 
find a convincing explanation to this trend. 

Besides the critical velocity, it is also essential to cal- 
culate the velocity reduction in terms of the system's pa- 
rameters. The velocity shift Av is defined as the differ- 
ence between the incident velocity and the velocity of the 
scattered soliton well after crossing the potential region. 
This depends on the depth and width of the potential well 
as well as on the magnitude of the incident velocity. Fig- 
ure[Hb) shows the shift in velocity as a function of the well 
depth Vo for fixed incident velocity v = 0.33. This veloc- 
ity is slightly larger than the critical velocity v cr = 0.3275 
for Vq = 4 and a = 2. The shift in velocity also depends 



on the magnitude of the incident velocity. The solid and 
dashed curves in Fig.UJa) show the shift in velocity for in- 
cident velocities greater than the critical velocity v cr and 
for fixed well depths Vo = 4.35 and Vq = 4, respectively, 
with a — 2. We calculate the velocity shift by monitoring 
the position of the density maxima as function of the evo- 
lution time. The slope of this curve gives the velocity and 
we calculate two velocities one before the interaction with 
the potential well and one after the interaction with the 
well. Both are measured far from the interaction region. 

To understand the physics of the diode's function, 
we employ a variational calculation that confirms our 
previous arguments based on the numerical results and 
specifically accounts for the velocity reduction shown in 
Figs. H (b) and H (a). Following Refs. 0113], we use an 
ansatz function that takes into account the incident soliton 
and a trapped mode, namely ^(x, t) = ij>i(x, t) + tp2{x, t), 
where ipi(x,t) = A\ sech[Ai (x — r\{)\ exp[i</>i + iv\x\ and 
ip 2 (x,t) = A 2 sech[x/qi] exp[i</> 2 + i a 2 \ogcosh(x/q 2 )]. 
Here A\, v%, rji, and </>i denote the amplitude, ve- 
locity, center of mass position, and the phase of the 
incoming soliton and A 2 , q 2 , a 2 , and 4> 2 denote the 
amplitude, width, chirp parameter, and the phase 
of the trapped mode, respectively. For simplicity, 
the potential is approximated by the delta function 
V(x) = — VqS(x). After constructing the Lagrangian, L = 

we derive the Euler-Lagrange equations for the variational 
parameters. From the set of resulting equations, we derive 
the following equation for the center of mass position of 
the soliton 

* = * - Jt (^) 771 " ) *' (4) 

where A\ /At = Vo A 2 sech( A\ r\\ ) sin(A0) and v\ = 
— Vo Ax scch(^4i 771) tanh(Ai 771) [A 2 cos(A0) + A\ scch(^4i r\\ 
with Acj) = <pi — <p 2 and ' denotes the time derivative. 
From this equation it is clear that the potential well plays 
the role of a damping factor with a damping coefficient 
—A\jA\. Assuming 2/3 be the effective distance over 
which the soliton experiences such a frictional force, the 
velocity reduction can be shown to take the form 

Aw ps -2V" /3AiA 2 sech(Ai r\ x ) sin(A0), (5) 

which is proportional to the trapped mode amplitude A 2 . 
If the transmission is perfect, i.e., A\ = 1 and A 2 = 0, 
there is no velocity reduction. From this simple argument 
it is obvious that the reduction in velocity arises primarily 
due to the imperfect transmission after the critical veloc- 
ity. In Fig. Qlb), the transmission coefficient curve ac- 
quires a slight downward curvature. It is this particular 
behavior that mainly accounts for the velocity reduction. 
Specifically, the coefficient of the trapped mode A 2 is pro- 
portional to 1 — T. To verify this, we have calculated A\ 
and A 2 numerically in terms of the incident velocity and 
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the well depth. In addition, we approximated A(j) by its 
value near the trap center, namely — .21tt (See Fig. 7 in 
Ref. [12), and set r\ = 0. Finally, we set Vb/3 = l/v 7 ^, 
which is a good approximation for delta potentials. As 
a result, we obtain the velocity reduction as a function 
of the well depth and the incident velocity plotted with 
dashed-dotted curves in Fig. EJb) and Fig. 0|a), respec- 
tively, which able to capture the main features qualita- 
tively. 

The mechanism by which the diode functions, can now 
be explained in Figure SJb). The solid and dashed curves 
show the reduced velocity of the soliton after crossing a 
single potential well for the two cases of Vq = 4.35 and 
Vq = 4 with a = 2. The dotted diagonal line shows the 
velocity of the incident soliton. The two horizontal lines 
mark the critical velocities of the two wells. For the left- 
to-right-moving soliton we follow the dashed curve marked 
as OABC ' . For OA, where the point A corresponds to 
v = v cr , the incident velocity is smaller than the critical 
velocity v l cr = 0.3275 of the left well and the soliton is 
reflected back to the left. For v > v l cr we follow the curve 
ABC, where the reduced velocity is well above the criti- 
cal velocity v r cr = 0.1175 of the right well and the soliton 
is transmitted to the right. Therefore, from left to right, 
the soliton is reflected for v < v l cr and transmitted for 
v > v cr . On the other hand, for the right-to- left- moving 
soliton we follow the solid curve marked as Oa(3^8. For 
Oa, the incident velocity is smaller than the critical ve- 
locity v' cr of the right well and the soliton is reflected back 
to the right. For v > v r cr we follow the line af3j, where 
the point 7 corresponds to v l cr + |Au r |. Here, the reduced 
velocity is well above the critical velocity v r cr of the right 
well but smaller than the critical velocity v\ r of the left 
well. As a result the soliton is reflected by the left well and 
transmitted to the right through the right well. For j6 the 
reduced velocity is greater than the critical velocity v l cr of 
the left well and the soliton is transmitted to the left. In 
conclusion, the right-moving soliton will transmit to the 
right for v > v cr and the left-moving soliton will reflect to 
the right for v < v l cr + \Av r \. As a result, the scattered 
soliton will always be flowing to the right irrespective of 
the direction of the incoming soliton provided its incident 
speed is in the velocity range v l cr < v < v l cr + |Au r |. The 
width of this velocity window \Av r \ ~ 0.01 agrees well 
with that obtained directly from the numerical solution of 
the GP equation, as can be seen by comparing Fig. [TJa) 
and the inset of Fig. EJb) . It is essential to note that the 
velocity reduction should be considerably larger than the 
velocity range over which the transition at the critical ve- 
locity takes place. If this condition is not met then partial 
transmission and reflection, with comparable magnitudes, 
takes place resulting in a failure of the function of the 
diode. 

When the two wells are well separated, the effect of the 
either of them is independent of the other. As we move 
the wells closer, they couple together after a certain small 
separation and start to affect the soliton propagation in a 



combined manner. Due to this coupling between the wells 
and nonlinearity of the soliton dynamics, there is a non- 
monotonous change in the velocity window for which the 
diode effect is observed. Figure[5ja) shows the width of the 
velocity window w as function of the separation s between 
the wells. We define the window width w as the difference 
between the critical velocities of the solitons transmitting 
from opposite directions. The separation s is the distance 
between the centers of the two wells. Figures [5jb) and 
(c) show the reflectance and transmittance as a function 
of the incident velocity of the soliton. Note that velocity 
window is now almost doubled compared to that shown in 
Figs.QJa) and (b). 

The directional flow of soliton with asymmetric poten- 
tial wells can also be achieved with two Gaussian wells 
with appropriate parameter combination. Gaussian po- 
tential supports both criticality |13j and velocity reduc- 
tion for soliton propagation. For example, the potential 
V(x) = -4.35exp[-3(a;- 6) 2 ] -4.0 cxp[-3(a; + 6) 2 ] shows 
diode behavior in the range 0.304 < v < 0.310. Therefore, 
the directional flow is not unique to a series of RM poten- 
tials. 

The combined effect of the velocity reduction and criti- 
cality can also be used to trap solitons. If the reduction in 
velocity changes the velocity such that it is smaller than 
any of the two critical velocities then the soliton will be 
trapped in between the potential wells and oscillates. One 
such parameters combination is v ~ 0.33, Vx = V2 = 4, 
a\ — 012 = 2, and x CTa \ = — ^cm2 = 6. For this set of 
parameters, if the soliton is incident from left or from 
right it can transmit through the first well and the re- 
duced velocity becomes v = 0.315, which is smaller than 
v cr = v r cr = 0.3275 and oscillates between the two wells. 

For multiple soliton cases, we find similar behaviors as 
for single soliton case only if the initial soliton separation 
is very large compared to soliton width. In this regard, 
we use a chain of two solitons, which are initially located 
at, say, x cm = 20 and x cm = 60. Qualitatively we get the 
same diode behavior as shown in Fig. [I] For smaller initial 
separation, the nonlinearity and the interaction near the 
potential regions destroy the diode effect. 

Discussion and Conclusions. Our numerical 
study of the dynamics of one-dimensional GP equation 
with two asymmetric potential wells (Gaussian or Rosen- 
Morse) shows that the soliton can propagate in one direc- 
tion for a certain velocity range and with certain param- 
eter combinations. This unidirectional propagation arises 
due to the criticality and the velocity reduction by a single 
potential well. When two such potentials with slightly dif- 
ferent parameters are in series, the combined effect causes 
the soliton to move in one direction only. A detailed ac- 
count of the critical velocities and velocity reductions of 
both potential wells explains this effect and predicts the 
velocity range over which the diode functions, in a good 
agreement with the numerical result. We also employ a 
variational calculation and estimate the velocity reduc- 
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tion as a function of the input velocity and the depth of 
the potential well. The semi analytical velocity reduction 
curves able to capture the main feature of the correspond- 
ing exact numerical ones qualitatively. 

The function of the soliton diode was demonstrated 
here for a single pair of potential depths V\ = 4.35 and 
V2 = 4.0. The width of the velocity window for which 
the diode works properly, certainly depends on the val- 
ues and ratio of these two parameters as well as on the 
widths of the potential wells and their separation. Thus, 
widening the velocity window and increasing the diode's 
efficiency requires further systematic investigation of the 
effect of these parameters as well investigating other types 
of potential wells. 

The GP equation considered here is the governing equa- 
tion of the dynamics of both matter- wave solitons in Bose- 
Einstcin condensates [2^1f3D] and optical solitons in optical 
fibers. Therefore, it is expected that the diode behavior 
found here to be realized in both cases. For the case of 
optical solitons, this diode may be a favorable instrument 
in achieving all-optical communication and computation 
in optical fibers. 
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Fig. 3: (Color Online) (a) Critical velocity a function 

of the well depth Vo of the RM potential with a = \/Vq. (b) 
Solid curve shows the velocity shift as a function of the well 
depth Vo of the RM potential with a = y/Vo for fixed incident 
velocity v = 0.33. Dash-dotted curve shows the velocity shift 
as a function of the well depth derived semi analytically. 



Fig. 1: (Color online) (a) Solid and dashed lines show the re- 
flection coefficients for the left and the right moving soliton, 
respectively, (b) Solid and dashed lines show the transmission 
coefficients for the left and the right moving soliton, respec- 
tively. The parameter set is as follows: Vi = 4.35, V2 = 4.0, 
«l = ai = 2, and rE cm i = — x cm % = 6. 
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Fig. 2: (Color online) Figures (a) and (b) show the spatio- 
temporal density plots with schematics for a case with \v\ = 
0.33, which corresponds to the dots in Figs. [Ha) and (b), re- 
spectively. The parameter set is same as in Fig. [T] 
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Fig. 4: (Color Online) (a) The velocity shift as a function of 
the incident velocity. Solid curve shows the velocity shift for 
Vo = 4.35 and a = 2 and dashed curve shows that for Vo = 4 
and a — 2. Dash-dotted curve shows the velocity shift as a 
function of the incident velocity derived semi analytically for 
Vo = 4. The vertical lines show the critical velocities for the 
two well depths, (b) Reduced velocity versus incident velocity. 
The horizontal lines show the critical velocities for the two well 
depths. The solid and dashed curves show the reduced velocity 
for the left and the right moving solitons, respectively. The 
inset shows, with the thick horizontal line, the velocity range 
for which the diode functions. 
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Fig. 5: (Color online) (a) Velocity window to as a function of 
the separation s between the centers of the wells. The width 
and depth of the wells are fixed. Dashed line is drawn to guide 
the eye. The points marked as A and B are the separations 
which are used for Fig. [T] and Figs. H](b) and (c), respectively, 
(b) Solid and dashed lines show the reflection coefficients for 
the left and the right moving solitons, respectively, (c) Solid 
and dashed lines show the transmission coefficients for the left 
and the right moving solitons, respectively. Position of the well 
centers are i cm i = —x crn 2 = 2. All the other parameters are 
same as in Fig. [T] 
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